A real quaternion spherical ensemble of random matrices 



One can identify a tripartite classification of random matrix ensembles into geomet- 
rical universality classes corresponding to the plane, the sphere and the anti-sphere. 
The plane is identified with Ginibre-type (iid) matrices and the anti-sphere with trun- 
cations of unitary matrices. This paper focusses on an ensemble corresponding to the 
sphere: matrices of the form Y = A X B, where A and B are independent N x N 
matrices with iid standard Gaussian real quaternion entries. By applying techniques 
similar to those used for the analogous complex and real spherical ensembles, the eigen- 
value jpdf and correlation functions are calculated. This completes the exploration of 
spherical matrices using the traditional Dyson indices /3 — 1, 2, 4. 

We find that the eigenvalue density (after stereographic projection onto the sphere) 
has a depletion of eigenvalues along a ring corresponding to the real axis, with re- 
flective symmetry about this ring. However, in the limit of large matrix dimension, 
this eigenvalue density approaches that of the corresponding complex ensemble, a den- 
sity which is uniform on the sphere. This result is in keeping with the spherical law 
(analogous to the circular law for iid matrices), which states that for matrices having 
the spherical structure Y = A _1 B, where A and B are independent, iid matrices the 
(stereographically projected) eigenvalue density tends to uniformity on the sphere. 

1 Introduction and main results 

This paper follows on from [22j where the author described a so-called geometrical tri- 
umvirate of random matrix ensembles; that is, classes of matrices that (in some sense) live 
naturally on the three surfaces of constant curvature: the plane, the sphere and the anti- 
sphere (this categorization was also used in [19] in a similar context). Through the circular 
law [HI [21 [T31 [32], we identify the Ginibre ensembles [12] (and other iid ensembles) with 
the plane, since in the large N limit the eigenvalue density is uniform in some planar disk 
and zero outside. In fact, in |31j it is shown that the general n-point correlation functions 
are universal in the class of iid matrices, assuming the first four moments are finite. 
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We identify the sphere with matrices Y = A B, where A and B are independent iid 
matrices. The spherical law (conjectured in and [29], and proven in [3]) establishes that 
the eigenvalues of these matrices (under stereographic projection) are uniformly distributed 
on the sphere in the limit of large matrix dimension. However, as yet nothing is known 
about the higher order correlations for the spherical ensembles beyond the Gaussian case. 
Lastly, the anti-sphere we may identify with truncated ensembles, since in the two studied 
cases (truncated unitary [331 123 M an d truncated orthogonal [18J) the eigenvalues are 
uniformly distributed in a hyperbolic space. There is, to date, no 'anti-spherical law'. 

This paper addresses the outstanding spherical ensemble in the traditional (3 = 1,2,4 
collection: the real quaternion (/3 = 4) spherical ensemble. Previously, the complex (/3 = 2) 
and real (j3 = 1) cases have been studied in [U HHJ EDI HSl E] • To quickly gain some insight 
into the behaviour of these eigenvalues we can simulate examples (see Figure [T]). Note the 




Figure 1: Plots of (stereographically projected) eigenvalues for 120 independent 100 x 100 
spherical matrices with real (left), complex (right) entries. 

ring of eigenvalues in the real spherical ensemble, which is absent in the complex ensemble; 
this ring corresponds to real eigenvalues, which have a finite probability of appearing in 
real ensembles. The remaining eigenvalues come in complex conjugate pairs, meaning that 
the real spherical ensemble has reflective symmetry across the ring (generally the complex 
eigenvalues of real matrices occur in complex conjugate pairs) and rotational symmetry 
about the axis through the centre of the ring. The complex ensemble has unbroken rota- 
tional symmetry. For the real quaternion spherical ensemble in the present study we will 
find that there is a depletion of eigenvalues along the ring corresponding to the real axis. 
Such an effect in the real quaternion Ginibre ensemble follows from the joint eigenvalue 
probability density function computed in Ginibre's original paper |12j . 

We note that the eigenvalue problem for matrices Y = A X B is equivalent to the 
generalized eigenvalue problem, that is solutions Xj to 

det(B - AA) = 0, 

where A is invertible. The real spherical ensemble was discussed in [7] in the context 
of these generalized eigenvalues, where the authors establish that the arctan of the real 
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eigenvalues are uniformly distributed on the unit circle. (See [TT] for more discussion on 
the statistics of the real eigenvalues.) 

To approach the real quaternion ensemble, we follow the skew-orthogonal polynomial 
method used in [11] , which includes an iterated integral technique (see Section [3]) used 
successfully for the (3 = 2 spherical ensemble [16] (a similar technique was also used in |12j). 
For a more detailed description of the techniques used in this paper and the underlying 
theory, see [22]. 

Since quaternions are crucial to this study, we first provide a quick overview. A quater- 
nion is analogous to a complex number, except that it has four basis elements instead of 
two. Typically they are written in the form q = qo + iq\ + jq<i + kq^, with the relations 
i 2 = j 2 = k 2 = ijk = —1, and the qi are in general complex. We will also use an alternative 
representation as 2 x 2 matrices: 



w x 

y z 



(1) 



where w = qo + iqi,x = qi + iqs,y = —qi + iq3,z = qo — iq\. The analogue of com- 
plex conjugation for quaternions we denote q* = qo — iqi — jq2 — kq^, or in the matrix 
representation 



— x 
w 



In the case that qo,qi,q2,Q3 G K we say that q is a real quaternion and from ([!]), with 
o = qo + iqi and /3 = q2 +iqs, we have 



a (3 
-P a 



with conjugate 



a 



n 



With Q = [qj^k] we denote by Q the matrix [q^j], and we call it the dual of Q. If 
Q = Q D then Q is said to be self-dual. We will regularly use quaternion analogues of the 
usual matrix trace and determinant [BJ. 

Definition 1.1. For eqi iVxJV matrix Q with real quaternion entries the quaternion trace 
is defined as the sum of the scalar parts of the diagonal entries 



N 

qTr Q := ^(<7o)j,j- 
i=i 



(3) 
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The quaternion determinant is defined by 



i 

qdet Q := ^ (-l) N ~ l J] qTr(q ab q bc ■ ■ ■ q sa ), (4) 
Pes N i 

where P is any permutation of (1,...,N) that consists of I disjoint cycles of the form 
(a — > b — > c —> • • • — > s — > a). 

Note that the definition ^ gives 

qTr QatxTV = ^Tr Q 2 iVx2Af, (5) 

where Q27Vx2JV is the matrix corresponding to QatxJV with the quaternions replaced by 
their 2x2 representatives ([2]). Furthermore, it is shown in [6] that with the definition Q 
and with Qat X 7V a self-dual real quaternion matrix 

* 1 /2 

qdet QatxJV = (detQ 2 Ar X 2Ar) • (6) 

Since we will be mostly using the 2x2 representation for the quaternions we will most 
often make use of ([5]) and ([6]) instead of Definition [lTJ although we will suppress the 'hat' 
notation for ease of reading. 

An element of the real quaternion spherical ensemble is a product of two real quaternion 
Ginibre matrices. These matrices, as defined in [12], are N x N matrices with independent 
Gaussian real quaternion entries. For ease of reference, we make an explicit definition. 

Definition 1.2. Define standard Gaussian random real quaternion variables as q = qo + 
iqi +312 + kQ3, where qo,qi, q2,Q3 are independent real standard Gaussian variables. 

Let A, B be independent N x N matrices with standard Gaussian real quaternion 
entries. Then the real quaternion spherical ensemble consists of N x N matrices of the 
form Y = A _1 B. 

Equivalently, choose A and B to be independent 2iV x 2N matrices containing the 
2x2 blocks Q, where a,/3 are complex numbers with real and imaginary parts being 
independent standard Gaussian variables. In this case Y = A _1 B is a 2N x 2N complex 
matrix. 

In Section [2] we calculate the matrix probability density function (pdf) for these ma- 
trices. It turns out that we can combine the results in [8] for the complex and real cases 
(also later derived independently in [TB] and [TT] respectively) to state the matrix pdfs for 
the three spherical ensembles (real, complex and real quaternion) in the one form. Note 
that in [8] the author shows that a matrix Cauchy distribution is obtained for all spherical 
matrices, independent of the elemental distribution. 
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Proposition 1.3. With Y a real = 1), complex = 2) or real quaternion = 4) 
spherical matrix, the matrix pdf is 

y_i r ( N+l+j g) 

Vp(Y) = vr-^ 2 / 2 TT \ 2 / det(l + YYt)-^, (7) 

where the 'dagger' should be interpreted as 'transpose', 'Hermitian conjugate' or 'quater- 
nion dual' respectively. Note that for (3 = 4, the determinant is to be understood as a 
quaternion determinant. 

We shall be interested in the eigenvalue distribution; in Figure[2]we present a simulation 
like those in Figure[TJ showing that (on stereographic projection) the eigenvalue distribution 
splits into two symmetrical hemispheres indicating repulsion from the ring corresponding 
to the real axis, with the complex eigenvalues coming in conjugate pairs (the latter is a 
general property of self-dual real quaternion matrices). 




Figure 2: (Stereographic) eigenvalue plot for 120 independent 100 x 100 spherical matrices 
with real quaternion entries. 

Using a Schur decomposition, in Section [3] we change variables in ([7]) to extract the 
eigenvalue joint probability density function (jpdf), by applying the iterated integration 
technique used in [16] and [TT]. Then, taking advantage of the rotational symmetry about 
one axis in Figure [2j we apply the fractional linear transformation 

A = (8) 

which takes the upper half plane to the unit disk D, and the lower half-plane to its com- 
plement, and we find the eigenvalue jpdf 
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with 



w = {wi, . . -,w N }, wj e 
1(1 11 1 



G 

W {Wi Wn J Wj 

r =( un(n-d/2 (*r N fr r(2iV + 2) 
1 J r(JV + i) jLl (r(2j)) 2 ' 

A(xi, ...,x n ) = ]~| (x k -Xj). 

l<j<k<n 

Note that the factors of (l/|x| — |x|) signify the expected repulsion from the real axis. 

The average over this jpdf is calculated in terms of a Pfaffian expression in Section 
[4j By finding appropriate skew-orthogonal polynomials (Section [5]), and using techniques 
similar to those in [3] we calculate the eigenvalue correlation functions in Section [6j finding 
the same structure as in the real quaternion Ginibre ensemble [12\ [T7] . 

In Section 6.1 we find the eigenvalue density to be 

1 (I _ r ) (2N + l)(2iV) ^ (2N - 1\ (r 2 ^- 2N - r 2 ^ 1 ^) 



rU r+ l) 2N + 2 vr J J 2N-2J-1 



P(l)( W > ~ ^2 , , ^2N+2 Z> 

3=0 

2 2N - 1 (2N + l)N f 1 -'' 2 



r 2 i i \ 2JV+2 - / 

r (r + i) 7T Jo 

where w = r e td , which for large iV tends to 



(l-T)"^ 1 /^!- 7/2)^-^7, (10) 



2N 

P(i)(w) ~ — — o\o ? (11) 

1 ; jv^oo 7r(l + r 2 ) 2 

which is the expected analogue of Eqn. (59)]. This limiting density tells us that, on 
stereographic projection, we have a uniform density of eigenvalues on the sphere in the 
limit of large matrix dimension. This is a consequence of the spherical law, which was 
discussed at the beginning of this introduction. 

With Xj i — ^ Xj/y2N in (22), then we note that (for large N) Q(X) approaches 

N 

C N He-^ 2 \X, -A/ H \X k - Aj| 2 |A fc -Aj| 2 , (12) 

j=l l<j<k<N 
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where Cm is independent of the Xj. This expression ( 12 ) is the same as the eigenvalue jpdf 
for the Ginibre ensemble of real quaternion matrices \12\ Eq. (2.28)], and so we may expect 
to recover the same bulk correlation functions (see [IT] ) from our correlation functions here 
by a suitable scaled limit. In Section 6.2 we apply the rescaling (39 ) and obtain Proposition 
6.2, which matches the result for the /3 = 4 Ginibre case. 



2 Matrix probability density function 

With A and B each being N x N matrices from the real quaternion Ginibre ensemble (that 
is, having iid standard normal real quaternion entries) the joint distribution of Y = A *B 
is 

(27r)- AN2 e-^ AAD+BBD y 2 (dA)(dB) = (2n)- 4N \-^ AAt+BBt » A (dA)(dB), (13) 

where we have changed to the 2x2 representation for the RHS. (For the remainder of 
this paper we will suppress the 'hat' notation, however we will endeavour to preserve the 
distinction between D for quaternion dual and f for Hermitian conjugate, which, along 
with 'qTr' and 'qdet', will serve to indicate which representation we are using.) The plan 
now is to change variables to Y and integrate out the remaining degrees of freedom. To 
achieve this we will follow the real quaternion analogue of the procedure used in [IB] and 
|llj . First we list some required lemmata specific to real quaternion matrices. 

Lemma 2.1 ([2B]). Let X be an n x m (n > m) matrix with real quaternion entries, and 
A = X D X. Then 

(dX) = c (qdet A) 2 ("- m ) +1 (dA), 

where c is independent o/X. 

Lemma 2.2 ([21]). With A and M both n x n real quaternion matrices (with A non- 
singular) and X = A D MA we have 

(dX) = (qdet A D A) 2n ~ 1 (dM). 

Lemma 2.3 QlOj). A corollary of the Selberg integral 130)1 is the evaluation 

poo poo N 

/ <&!••• / dx N ]Jx 2a e- 2xi J] \ x l~ x it 
"'° ^° 1=1 i<i<i<N 

N-l 

= 2~ 2N ( N+a 1 []r(3 + 2j) T(2a + 1 + 2j). 

3=0 
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Lemma 2.4. Let B mxn = x mxm A mxn y nxn , where each of the matrices consist of 2 x 2 
real quaternions then 

(dB mxn ) = |qdet x mxm | 4n |qdet y n xn\ 4m (dA mxn ). 

This last lemma is the real quaternion analogue of the statement about real matrices 
in Theorem 2.1.5 of |24j . From the working there we can verify this lemma by noting that 
x acts on the columns of A (of which there are n) and y on the rows (of which there are 
m), and then each quaternion has 4 independent components. 

Letting C := AA D (which is positive definite) we see from Lemma 



2.1 



that 



(dA) = cqdet C (dC), 



(14) 



and also, letting B = AY, by Lemma 2.4 we have (a!B) = |qdet A| 4Ar (<iY) 
= (det C) N (dY), and so (13) becomes 

(27r)- 47v2 e- Tr ( c ( 1+YYt ))(detC) iV+1 / 2 ((iC)( ( iY). 



(15) 



To calculate c we multiply both sides of (14) by e ( TrAAt )/ 4 a nd integrate over the whole 
domain 



- (TrAAt )/4 



(dA) 



-(TrC)/4 



C>0 



(detC) 1/2 (dC), 



where the positive definiteness of C is indicated by C > in the integral terminal. We 
note that the LHS is just an integral over AN 2 independent standard Gaussians and so 



(27T) 



2N 2 



/ c>0 e-M/4(detC)V^C)' 



Let D := (1 + YY t ) 1 / 2 C(l + YY+) 1 / 2 in ^ and then integrate over D to obtain 

♦ ™rfn^e- ( ™ )/4 (detD)^/2(dD) 

(2ir)~ 2N det 1 4- YY"h~ 2jV D>0 - (dY) 

[ ] d6ti + J J c>0 e-(T r c)/4 (detC )i/2( d c) {d h 



where we used Lemma 2.2 to obtain the Jacobian for changing variables from C to D. 
Since both C and D are self-dual real quaternion matrices we can make use of the well 
known Jacobians for changing variables from the matrix entries to the matrix eigenvalues 
(see for example |10| Chapter 1.3]) so that the ratio of integrals becomes 

J™ dXr... J™ d\ N nf =1 Af + W 2 n j<k \\ k - Aj f 



J Q oo dx 1 ...j oo dx N uU x ^- Xj/2 H. 



j<k 



|Afe — A ? - 1 4 



S 



By changing variables Xj t— >■ 4Xj we can apply Lemma 2.3 to evaluate the integrals, giving 
the final expression 



7T 



-2»^ 1 r(2(JV + i+j)) 



n 

3=0 



r(2(i + i)) 



qdet(l + YatY 



By comparing this to the analogous results in the real and complex cases we can write 
them in the single equation ([7]). 



3 Eigenvalue jpdf 



Here we will change variables in the matrix pdf from the elements of Y to the eigenvalues 
of Y. Using the 2x2 representation of quaternions Q apply the Schur decomposition 
Y = Q^R^atQ (see [JJ for a detailed description of this in the real quaternion case), with 
R-2JV an upper triangular real quaternion matrix 



R 



■2iV 



Ai 

Ai 







ri,2 Sl,2 

_0 A 2 
A 2 



-si,n Fi,jv 

T2,N S 2 ,N 
~S2,N r 2 N 



" 




"00" 



















Aiv 






where Xj and Xj are the N complex conjugate-paired eigenvalues of Y, Tj & and Sjk are 
complex numbers, and Q E Sp{N) / (U {1)) N is a symplectic (ie. unitary real quaternion) 
matrix (using the 2x2 representation). Since the decomposition is only unique up to 
ordering on the eigenvalues, we divide through by 2 N T[N + 1]. The Jacobian for this 
change of variables is [231 E] 

(dY) = (QtdQ)(dR 2JV )|A(A)| 2 , 

where A(A) := fX^li^j — ^j) Y\i<j<k<N(^j ~ ^k)(^j ~ ^k) and (dX) represents the wedge 
product of all the independent elements of a matrix of differentials dX. Then 



V 4 (Y)(dY) 



N-l 



(QtdQ)(dR 2jV ) " - r(2(JV + l+j)) 
„2N*2 N T(N + 1) f} Q r (2(j + 1)) 



|A(A)| 2 det(l 2 7v + R27vR 2W ) 



-2N 



(16) 



We will now apply the iterated integration technique used in [16] and [TT] to integrate 
over the strictly upper triangular components of R2N- First we separate out the last 
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column 



R 



■27V 



R2JV-2 U2N-2 

z N 



where z^ is the bottom right 2x2 block on the diagonal, which contains the eigenvalues 
Aat, Ajv, and U2N-2 is of size (2N — 2) x 2; explicitly 



U-2N-2 



ri,N-l 

-Sl,N 



Sl,N 



rN-l,N-l Sn-1,N 
— SN-1,N rN-l,N 



Then the determinant in ( 16 ) becomes 



det(l2jv + R2JVR2JV) 2N ~ det 



12AT-2 + R2JV-2R.2JV-2 + U2N-2U 2N _ 2 U'2N-2Z N 

1 2 + ZpfZlf 



det 



Z N u 2N-2 
t 



+ R 2 JV-2R2iv_2 + U 2N~2U\ N _ 2 — Uz' N (l + Z N z' N ) Z N u\ N 
z NU 2n _2 





1 2 + ZNZ 



= det (l 2 + ZNZff) det ^l2JV-2 + R2JV-2R2]V-2 + u 2N-2u\ N _ 2 — U2jv(l + ZnZ^) 1 Z N u\ n _,^ 

We use the identity z*(l + zz*)z = 1 — (1 + z^z)^ 1 to write the second determinant above 



as 



det U2JV-2 + R27V-2R2W-2 + u 2/v-2(l2 + ZnZn) 1- 4tv-2 

= det ( l 2 JV-2 + R27V-2RLV-2) det + v 2N-2 V 2N-2) , 



where 



V2N-2 ■= (12JV-2 + R2Ar-2RLv_ 2 ) ^ ^2iV-2 (12 + Z N Zn) 1/2 . 



(17) 



So (16) becomes 



p 4 (Y)(dY) - n2N * 2Nr{N + 1) 11 r 



2.V) Vr(2(iY + l+j)),x 



! I) r(2(i + i)) 



|A(A)| 2 det(l 2 + Z7v4 



-2N 



X det ( 12AT-2 + R2Af-2R2AT_ 2 

Note that 



-2N 



det ( 1 2 + v! 2N _ 2 v 2 n-2 



-2N 



(18) 



AT AT-1 
(dR 2A r) = Y[ dXj f\ (du 2N-2s) 

3=1 s=l 
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and we want to integrate over the independent elements of U2N-2, ie. we want to calculate 

-2AT 



J (du 2 N-2) det (l 2 + v\ N _ 2 V2N-2 



(19) 



Applying Lemma |2.4| to (17) we have 

(du 2N -2) = qdet(ljv_i + R JV _iR^_ 1 ) 2 qdet(l + z N z%) 2N ~ 2 (dv 2N -2) 

= det(l 2 AT-2 + R27V-2R-2AT-2) det ( 1 2 + Z N Z ] N ) N ^ {dv 2N - 2 ) 



(20) 



For ease of notation, let CL2N-2 '■= v\ N _ 2 v< iN-2, so 



0-2N-2 



s^N-l I |2 1 

L 7 =i \ r iM +\sj, N \ 









22j=l \rj,N\ + 



We can again make use of Lemma [2TT] by applying it to «2Af-2 giving 

(dv 2 N-2) = c (qdet a 2 7v_2) 2iV ~ 3 (da2Ar-2) = c (det a 2 iY_ 2 ) Ar " 3/2 (da2Ar-2) 



(21) 



Multiplying both sides of (21) by e Tr ( v 2N-2 v ™ 2) anc j integrating we have for the LHS 



-Tr{vl N _ 2 v 2 N-2) 



(dv 2 N- 



I, 



2EjL - i 1 (l«i,Ar[ 2 +|r J) Ar| a ) 



iV-1 



(-oo,oo) 4JV ~ 4 
4iV-4 

- n 

and for the RHS 



"fj dRe(sj i Ar)(iIm(sj i Ar)dRe(rj j 7v)rfIm(rj i Ar) 



e 2x2 



7j- \ 2N—2 
2/ 



c / e- Tra ^- 2 (deta 2 7V-2 



,iV-3/2 



(do,2N- 



X 2x^AT-3^_ s r(2JV-2) 



4JV-1 



because CL2N-2 = diag(x,x), where x G (0,oo). We deduce that 

^2Af-2 



7T" 



r(2iV - 2)' 
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/ 



Now to calculate (19) 

I t N ~ 2A 

(du 2 N-2) det ( 1 2 + v ] 2N _ 2 v 2 n~2 



det(l 2 AT-2 + R 2 Ar-2R.2Ar_ 2 ) det(l 2 + znz^) 1 *' 1 I (dv 2 N-2) det ( 1 2 + v\ n _ 2 v 2 n-2 i 
K$&-2 T ?2N- 2) I det ( a 2N~2) N ~ 3/2 det(l + a 2N . 2 )- 2N (da 2N _ 2 ) 



JX 2N-2 I ■' ' ax — IX 2N-2 I \l + x J (l _|_ x yN+3 " X ' 







27V-3 



1 



where -^ 2 jv_ 2 stands for the determinants in front of the integral after the first equality, 

and i4S-2 = ^2^-2 x vr 27V - 2 /r(2iV - 2). Letting y = x/(l + x), and so dx = (1 + x) 2 c?y, 
we find a Selberg integral |30J resulting in 



\ 2iV-3 -t r l 

/V 2.\-:> / + n + ^2^+3 ° X ~~ ^2iV-2 / ■' / M ' ' "'^ 







(2) r(2iV - 2)r(2iV + 2) 
2iV - 2 T(AN) 



7r2Ar _ 2 r(2|v + 2) det(l2iv _ 2 + R2Ar _ 2R t 2jv _ 2 ) det(l2 + 



Substituting this into (18) we have 

: v mrjvl _ ^- 2(JV(jV - 1)+1) r(2iv + 2) ^ r (2(iv + 1 + j)) | A(X)| 2 
1 A J 2 Jv r(iv + i) r(4iv) 1JL r(2(j + i)) (i + ia^i 2 ) 2 ^ 1 ) 

X det (l 2N - 2 + R 2 Ar- 2 R2iV-2) II ^ A ( du 2N-2s) (Q^Q), 

j=l s=2 

where the * on the integral indicates that we have only done some of the integrals required. 
We then iterate this procedure over each of the N — 2 remaining columns. 

Lastly, we need the integral over (QtdQ). As mentioned above, Q G Sp(N)/(U(l)) N , 
and using a result from [25] (see \22\ Ch. 7.1] for more details on this point) we have 



t _ V1W2AQ) _ _^ 



(Voi(f;(i)))« nf,ir(2j)' 

where we have used the result from Eq. (108)] 



Vol (Sp(N) 



2 N n N(N+l) 
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So having integrated out all independent variables in (18), except for those correspond- 
ing to the eigenvalues Xj = Xj + iyj, we are left with the eigenvalue jpdf 

7T~ N ^ T T{2N + 2) 



N IA - A I 2 N 

x n ?TTtVpwi) n iA fe -A/iA fc -A/n^<%> ( 22 ) 

j = l v ' ' l<j<k<N j=l 

where, for later convenience, we restrict A to the upper half plane (in other words we order 
the complex conjugate pairs Im(A) > Im(A)) multiplying through by 2^. 

3.1 Fractional linear transformation 

As in [H] we can take advantage of the expected rotational symmetry by applying the 
fractional linear transformation Q. Then with Aj = Xj + iyj and Wj = aj + ibj we have 

Vj = -. ■ — ^-L and dxjdyj = -. r-rdajdbj, 

1 1 + W j | 1 |l + U>j| 4 

in which case 



n iA fc -A/iA fe -A/= n \ w . _i_ ii^L + - w j\ 2 \ w J w k - m 2 

l<j<k<N l<j<k<N 1 3 



Also 

2 



Substituting into ( 22 ) and simplifying we have 



e(w)(*r) = ( -i)^w^ n 



ir)- N -£- r(2A^ + 2) 

^+i)l = l (r(2j)) 2 

N 1 /in \ N—l/2 1 -lio-l / ,\ 



from which ([9]) follows, where A(x) is the standard Vandermonde product. (Note that the 
factor of (i) in Q is an artifact of splitting the product l/|iOj| — \ wj\ into the factors of 
r(wj) and r(l/wj). It was for this splitting that we ordered the complex conjugate pairs 
in (ph.) 
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4 Generalized partition function 

We will now find a Pfaffian expression for the generalized partition function 



N . 
Z n[v] := (Yl v i w j)) Q = J V ( Wl ") ( dwi ^ 



v(w N ) (dw N ) Q(w), 



(24) 



where O is the unit disk, with which we shall calculate the correlation functions in Section 
[6} Using the Vandermonde identity we have 



A w, 



w 



det 







j=l,...,N 






j=l,...,N _ 



k=l,...,2N 



= { _ 1) N(N-l)/2 det 

and so ^ becomes 

Q^) = (-l) N{N - 1)/2 C N {l-^— 2 r( 



Pfc-l(=r) 



i=i,...,N 

fc=l,...,2JV 



[Wj)T 



3=1 



det 



Pk-i{=.) 



(25) 



j=l,...,N 



Before proceeding we will reorder the columns of the determinant in (25), as in [IT, eqn. 
(30)], as follows 

P0,PN-l,P2,PN-3, ■ ■ -,PN-2,Pl, (26) 

introducing another factor of (— 1) W 2 )^/ 2-1 )/ 2 . For labeling purposes we define 

Q2j(x) = P2j(x), Q2j+l(x) = p N -l-2j(x). 

This reordering will make finding the relevant skew-orthogonal polynomials straightforward 
(see Section [5]) . With the relabeling we have 



N 



Q(W) = C N X\-^T( Wj )T 1 



3=1 



W , 



det 



qk-i(wj) 
gfe-l(=r) 



(27) 



3=1,-- ..JV 
k—l,....2N 



Substituting this into (24), expanding out the determinant using the Leibniz formula, 



and then ordering the terms to give a Pfaffian we obtain the following result (see [TO 
Chapter 6.1.2] for explicit details). 
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Proposition 4.1. The generalized partition function for the real quaternion spherical en- 
semble can be written as 

Z N [v] = T[N + 1] C N Pf[ lj>k [v]] j;k=lj ,.. t2N , (28) 

where 

ljM := \ L w T{w)T \h) (k) " qj ~ i (k) qk - i{w) ) {dw) - 

5 Skew-orthogonal polynomials 

In this section we find the monic polynomials that will reduce the Pfaffian in ( |28| ) (with 
v = 1) to block diagonal form, where the diagonal blocks are the 2x2 matrices 

hj ~ 


in which case 

N 

Z N [l]=r[N + l]C N ]Jh j . (29) 
i=i 

In other words, define the skew-symmetric inner product 

(q.j,qk) ■= 7j+i,fe+i[i] 

and look for polynomials to satisfy the skew-orthogonality conditions 

(P2j,P2k) = (P2j+l,P2k+l) = 0, (P2j,P2k+l) = ~(P2k+l,P2j) = Sj,k h j- (30) 

By writing out 7_j fe[l] explicitly, we see that the only time it is non-zero is when the 



angular dependence cancels. With this fact, and the convenient reordering in (26), it can 
be seen that the skew-orthogonal polynomials are the monomials 

q 2j ( w ) = q 2j +l(w) = w 2N - l - 2 K (31) 



With w = re tB we substitute these polynomials into 7j,fc[l] to obtain 

"1 l_ r 2 

o (1 + ^ 2 

By changing variables c = r 2 we can make use of the identity in [15, Equation 3.216 (1)]: 
for general a, b such that Re b > 0, Re (a — b) > 0, 



hi = 1*1 n J,2 N+ 2 ^ l -r iN -^)dr. 



f (t b ~ 1 + t a ~ b - 1 )(l + t)- a dt 
Jo 



r(&)r(a - b) 

—YioT) 
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and so, with b = y, a — b = x, we have a non-standard form of the beta integral 



/V>(i-*r'*=S^. 

Jo r(x + !/) 



This gives 



vr(2iV_4j - 1) /27V- IV 1 



and, from (28), 



(2JV+1)(2JV) V 2j 



JV-1 
j=0 



cwro/v + i)" 



j = 0,...,N-l, 



(32) 



6 Correlation functions 



To calculate the n-point correlation functions we will use the functional differentiation 
formula 



P(n)(w) 



1 



6" 



Zn[v] 5v(w\) ■ ■ ■ 5v(w r 



-Z N [v] 



v=l 



(33) 



Indeed, by the skew-orthogonality (30), it is straightforward to calculate the 1-pt correlation 



function (the eigenvalue density) by applying (33) to (28) 

N ( d 



Zn[1] 



3=1 



h) n h h 



fc=0 

k^j-1 



which reduces to (10) by making use of (29) and (32). 



For the general correlations, we make use of the formula det(l+AB) = det(l + BA), or 
its quaternion determinant analogue qdet(l + AB) = qdet(l + BA) (for BA self-dual) - 
or an equivalent Pfafhan formula, see [28J — in conjunction with the Fredholm quaternion 
determinant 

00 roo poo 

qdetfl + XK] := 1 + ^ — / dx 1 ■ ■ ■ dx s qdet[K(xj,x k )]j,k=i,...,s, 

°" J — CO J — oo 

where K is an integral operator with 2x2 kernel K(x,y), and the matrix [K(xj , j,fe=l,...,s 
is self-dual. The details are lengthy and tedious, however it is directly analogous to the 
same approach for the real spherical and real Ginibre ensembles, see [22] for a detailed 
exposition (the techniques therein were inspired by those in [1] and |10j). We find that the 
correlation functions are of the same structure as those for the real spherical ensemble [11] . 
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Proposition 6.1. With the eigenvalues transformed according to Q), then using the skew- 
orthogonal polynomials (31) the n-point eigenvalue correlation functions for the real quater- 
nion spherical ensemble are 



P( n )(wi, ...,W n ) = Pf [K N (wi,W m )] , Wi £ 



K N (x,y) 



D(x,y) S(x,y) 
-S(y,x) I(x,y) 



(34) 



(35) 



where 



N-l 



D(x,y) = ^2 ^{ a 2j(x)a 2j+ i(y) - a 2j+ i( 



j=0 



S{x,y) = - ^2 j-,(a2j{x)b 2 j+i{y) - a 2j+1 (x)b 2j (yfj , 

j=0 3 

1 JV ~ 1 1 

I(x,y) = --^2 J^{ b 2j( x ) b 2j+i{y) - b 2j+1 (x)b 2j (y) 



=o "3 



and 



aj(x) = \x\ 1 r(x)qj(x), 
b j (x) = \x\- 1 r^ 

Substituting in the various factors we have 

(2N + 1)(2N) ^ (2N - 1\ {(wzf3+V 2 - N - (u,*)*" 1 /*-^ 



S(w, z) = Ay 



7T 



E 

j=Q 



2N - 4j - 1 



An 



(2N + 1)(2N) ^ (2N - 1\ ((wz)3+y 2 - N - {wz) N - 1 / 2 ^) 



7T 



N-l 

E 

3=0 



where 



1 



1/2 



2N - 2j - 1 



1/2 



(36) 



7V+1 



r z + 



N+l ■ 
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with I 

w| — Tin and \z\ — t z . The other kernel elements satisfy the relations 

I{w,z) 
D(w,z) 



\w 

1 



w 



S w 



1 



(37) 



In the following, we analyze the eigenvalue density, which includes several checks of our 
working, however we first rewrite S(w,z). Let 



JV-l 

cr( 7 ) := £ 
3=0 



2iV — 2 j" — 1 



then 



7- 



,dff(7) ^ (1 + 7) 
c?7 



2JV-1 



2 7 AT-l/2 ' 



by anti-differentiating we have 
a(l)-a(v)-- 



l da('T) J 1 /' 1 (l + 7) 



\2N-1 



drf ~' 2 J v 7 n+i/2 
From the definition of £7(7) we see that <r(l) = and therefore 



-dry. 



a{v) = 2 



2N-2 



1-v 



(l-lT N - 1/2 (l- 7/2)^7, 



where we have changed variables 7 1— )• 1 — 7, thus giving the integral expression 
S(w, z) = A j 2N ' l{2N + l)N r W \l - 7 )-^(l - 7/2)^7, 

7T JO 



(38) 



which we will make use of in analyzing the asymptotic behaviour of the correlation func- 
tions. 



6.1 Eigenvalue density 



With n = 1 (34) reduces to /)m(io) = S(w,w). We see from (36) and (38) that this 
immediately gives us the two expressions for the density in (10). We can check that 



2tt pi 

d6 I rdr P{1) (re ie ) = N 
Jo 
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Figure 3: The points are the bin mid-points for a simulated radial (normalized) density 
of eigenvalues of 2000 instances of 500 x 500 matrices, where the eigenvalues have been 
transformed according to Q. The noisiness of the data near r = is due to the relatively 
small number of eigenvalues near the origin, which is a consequence of the relatively smaller 



area one is sampling from. The solid line is a plot of (10), scaled by 1/N. 



as required. For another check, in Figure 3 we find good agreement between ( 10 ) and 



a simulated eigenvalue density. Note the elevated density near r = 1, which we recall 
corresponds to the real line. We interpret this as the effect of repulsion from the real line, 
as discussed in the Introduction. 

To analyze the behaviour for large N we first exponentiate the logarithm of the inte- 



grand in (10) giving 

»l-r 2 



I 

JO 



exp [(2N - 1) log(l - 7/2) - (N + 1/2) log(l - 7)] dry. 



For large N this integral will be dominated by points close to the upper terminal, so we 
Taylor expand the exponent about 1 — r 2 



(2N - 1) log(l - 7/2) ~(N + 1/2) log(l - 7) 

( 1 + r 2)27V-l\ ^ ^ [2N+1 



log 



227V-l r 27V+l 



+ (7 " 1 + 



2r 2 



2N- 1 
1 + r 2 



and so for large N the integral in ( 10 ) is approximated by 

„2\2AT-1 rl-r 2 



(1 + 7^ 



22AT-l r 2A r +l 



exp 







(7 - 1 + r l 



2iV + l 2N -1 



2r 2 



1 + r 2 



d'y. 



Substituting this into (10) and taking N to infinity (noting that r — 1 < 0) yields (11). 
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6.2 Scaled limit 



Here we aim to recover the bulk real quaternion Ginibre correlation functions of [T7], by 
zooming in on the origin. Under the transformation ^ the origin in the original coordinates 
is mapped to the point (1,0) and so we make the replacements 



w h> w * = 1 + 2iW 
z ^ z tt _ i + 2iZ 



7T 



W = X + iY, 
Z = A + iB, 



(39) 



in ( 38 ) and look for the large N limit of 



47T 

iV 



S(w*,l&) = 2 27V+1 (2iV + l)A 



(l_ 7 )-^-i/2 (1 _ 7/2 )2iV-i d7 . (4 ) 



We have chosen the transformations ( 39 ) to make a correspondence with the results in |17| , 
however we note that the density will be twice as large as in that work since here we have 
restricted the eigenvalues to be in the upper half plane. 
First rewrite A,,,t v i as 



AT. 



(-1 _ r 2 \l/2/i _ 2 \l/2 
-1/2 I 1 r »;^ V 1 ^ 

(l + ^+Hl + ri)^ 1 ' 



where 



Using the asymptotic behaviour 



, AJirY A \WV 
AT iV 



2\ 1/2 



1 + 



xVN-x 2 /2 



~ e 

which holds up to terms of order 1/y/N in the exponent, we find that 

r^~ 1/2 ~ exp f-2yv / vr]V + 2vr|H/| 2 - 4vry 2 ) 

2vr|W r | 2 - 2VttNY- 2nY 2 



(41) 



(l + r^+^^exp 



(1 -r 



2 .W2 



i/2 f^y 



N 



1/2 
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and so 



7T (YB) 



1/2 



-27r(y 2 +B 2 ) 



Ignoring terms of order 1 /iV we take 



and the integral in (40) becomes 



I 



(1-7)' 



JV-1/2. 



1-7/2 



,2JV-1 



c?7 



1 + 



7 



2\ N 



since, for large N, 7 is restricted to lie close to zero. Then, by changing variables 7 1— > 
—2ij/\/N the integral approaches 



7T 



erf (y^W - Z)) . 



Combining the previous we have the large N behaviour of (40). 



Proposition 6.2. Changing to the variables w* and z* from (39) we find (38) has the 
limiting behaviour 



%S{w\z*) ~ -(YB)^ 2 e- 2 ^ Y2+B ^erf(^(W-Z) 
which, with W = Z , gives us the scaled limit of the density 

e^ y2u2 du. 



(42) 







IGttYV 4 ^ 

The scaled limit of the remaining kernel elements (D and I) can be obtained similarly, 
or from (42) by using formulae analogous to (37) that apply to the real quaternion Ginibre 
ensemble 

I(w, z) = iS(w, z), 
D(w,z) = —iS(w,z). 



The result from [T7j corresponding to ( 42 ) is 
2vr 



(YB) l ' 2 e-^ w \ + \ z \^ z +w htf(^{W-Z)). 



(43) 



Ignoring the imaginary parts of (43), which always cancel in Kanzieper's formulation, we 
can use the properties of the complex error function to see that the correlation functions 
for the real quaternion spherical ensemble agree with those of the corresponding Ginibre 
ensemble under the scaling limit implied by (39) (up to the aforementioned factor of 2). 
Further, from the recent results in [31] . we see that in this limit, we recover the universal 
correlation functions for planar (iid, with first four moments finite) ensembles. 
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7 Discussion 



As discussed in the introduction, a tripartite categorization of matrix ensembles into 
geometrically-themed universality classes was suggested in [22J (which was also used in 
a similar context in [E]). Each of these classes corresponds to a surface of uniform curva- 
ture: the plane, the sphere and the anti- or pseudo-sphere. The Ginibre ensembles, being 
ensembles of iid matrices, obey the circular law, which states that the eigenvalues are uni- 
formly distributed on a planar disk in the limit of large matrix dimension, and so we may 
associate the Ginibre ensembles with the plane. 

Taking two of these iid matrices and forming the product Y = A _1 B, which is a matrix 
generalization of a Cauchy random variable, gives us spherical ensembles. The spherical 
law (analogously to the circular law) states that the eigenvalues have uniform density on 
a sphere (under stereographic projection) in the limit of large matrix dimension and so 
these ensembles are naturally associated with the sphere. In this work we have applied the 
(skew-) orthogonal polynomial method to calculate the correlation functions for the real 
quaternion spherical ensemble, which is the ft = 4 analogue of the ensembles considered 
in |X9(, [20l [T6] and [11]. This completes the exploration of the matrices in the spherical 
universality class using the Dyson indices ft = 1, 2, 4. As the spherical law was established 
in [3] by adapting the proof of the circular law in [32] , we suspect it is possible to adapt 
the proof in [31J to show that the general n-point correlation functions for the complex 
spherical ensemble are universal for the spherical class of matrices. 

The third universality class in this categorization corresponds to the anti-sphere, which 
is a surface of constant negative curvature. The ensembles in this class are given by 
sub-blocks (or truncations) of unitary matrices. To date, truncations of complex unitary 
[331 GH E] and real orthogonal p2] matrices have been analyzed (for a discussion of the 
real orthogonal case also see [22]). Naturally enough, the next step would seem to be an 
analysis of a ft = 4 real quaternion analogue of truncated unitary and orthogonal ensembles. 
Further, to conclude the program laid down in [22] . perhaps by again adapting the proof 
of the circular law in [32], then a proof of a so-called 'anti-spherical' law may be achieved 
(see [22, Chapter 7.6.2] for more on this), and, even further, perhaps using the methods of 
|31] find that the general correlations are universal. 
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